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Recent observations suggest that many old pulsar wind nebulae (PWNe) are bright TeV 7-ray 
sources without a strong X-ray counterpart. In this paper, we study the spectral evolution of PWNe 
taking into account the energy which was injected when they were young for old PWNe. We model 
the evolution of the magnetic field and solve for the particle distribution inside a uniformly expanding 
PWN. The model is calibrated by fitting the calculated spectrum to the observations of the Crab 
Nebula at an age of a thousand years. We find that only a small fraction of the injected energy from 
the Crab Pulsar goes to the magnetic field, consistent with previous studies. The spectral evolution 
model of the Crab Nebula shows that the flux ratio of TeV 7-rays to X-rays increases with time, 
which implies that old PWNe are faint at X-rays, but not at TeV 7-rays. The increase of this ratio 
is primarily because the magnetic field decreases with time and is not because the X-ray emitting 
particles are cooled more rapidly than the TeV 7-ray emitting particles. Our spectral evolution 
model matches the observed rate of the radio flux decrease of the Crab Nebula. 



I. INTRODUCTION 

A pulsar releases its rotational energy as a relativis- 
tic magnetized outflow called a pulsar wind. The pul- 
sar wind collides with surrounding supernova ejecta, 
forms the termination shock, and creates a PWN [9(. 
The acceleration of the pulsar wind particles occurs 
at the termination shock and the PWN consists of 
the magnetic field and the ultrarelativistic particles 
[Tol Il2j . Such a created PWN emits photons ranging 
from radio to TeV 7-rays via the synchrotron radi- 
ation and the inverse Compton scattering. Current 
status of theoretical models as well as observational 
confrontations is reviewed by Gaensler & Slane (2006) 

a- 

The Crab Nebula is one of the best studied PWN 
at almost all observable wavelengths including its cen- 
tral pulsar, called the Crab Pulsar. Kennel & Coroniti 
(1984) § studied the spatial structure of the Crab 
Nebula, assuming that it is a steady state object (KC 
model). They found that the magnetization parame- 
ter a, the ratio of the electromagnetic energy flux to 
the particle energy flux just upstream the termination 
shock of the pulsar wind, must be as small as 0.003 to 
explain the observed dynamical properties of the Crab 
Nebula. Atoyan & Aharonian (1996) [l| succeeded to 
reconstruct the current observed broadband spectrum 
of the Crab Nebula by the use of the KC model. 

The KC model is not fit to consider the evolution 
because it is a steady state model. However, it is very 
important to consider the spectral evolution of the 
Crab Nebula. We need to consider the spectral evolu- 
tion to explain the flux decrease rate of the Crab Neb- 
ula at radio and optical wavelengths 0, EH, EH . More- 
over, we need to understand the spectral evolution 
of PWNe. Recent observations found many PWNe 
which have a variety of the characteristics, such as the 
ages, the expansion velocities, the morphologies, and 
the spectra of PWNe. This variety must be related 



to the spectral evolution of PWNe. For example, de 
Jager & Djannati-Ata'f (2008) @ discussed the pos- 
sibility that the aged PWNe can be seen as the TeV 
7-ray sources without a strong X-ray counterpart. 

In this paper, we revisit a spectral evolution model 
of PWNe. Although several spectral evolution mod- 
els of PWNe have been studied d, 0, E3] , several is- 
sues still remain to be further clarified. In our model, 
the PWN is simply treated as an expanding uniform 
sphere. We do not take into account effects of the 
spatial structure of the PWN, since it is somewhat 
costly and too detailed to see the whole spectrum. 
The energy inside the PWN is injected from the pul- 
sar spin-down energy, which is distributed between 
the relativistic particles and the magnetic field with 
a constant ratio. Our simple model can describe the 
observed basic features. 

We study the spectral evolution of the Crab Neb- 
ula for the first application of our model. The Crab 
Nebula can be used as a calibrator of our model when 
it will be applied to other PWNe in future. In section 
HT1 we describe our model of the PWN evolution. In 
section IIII1 we apply this model to the Crab Nebula. 
In section HVl discussions and conclusions are made. 



II. THE MODEL 

A. Energy Injection 

In this paper, we assume that the PWN is a uniform 
sphere expanding at a constant velocity upwn- The 
assumption of a constant velocity i>pwn is the easiest 
way to take into account the expansion of the PWN, 
although the real behavior of the expansion velocity 
could be more complex as the work made by Gelfand 
et al. (2009) [7]. We consider that the age of the 
PWN is younger than lOkyr in this paper and that the 
constant velocity i>pwn would be a good assumption 
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in this range. 

For the components inside the PWN, we assume 
that the PWN is composed of only magnetic field and 
rclativistic electron-positron plasma, both of which 
are injected from the pulsar inside the PWN. We 
divide the energy injection from the pulsar into the 
magnetic field energy E^ and the relativistic particle 
energy E c using the time independent parameter rj 
(0 < i] < 1). The fraction parameter n is the injection 
ratio of the magnetic field energy to the spin-down en- 
ergy. The fraction parameter rj in our model is similar 
to the magnetization parameter a in the KC model, 
although they are not the same. The magnetization 
parameter a is the ratio E^/E e at the pulsar wind 
region immediately upstream the termination shock. 
On the other hand, the fraction parameter rj pertains 
to Eb/(Eb + E c ) into the PWN region uniformly. 

For the particle injection, we also need to deter- 
mine the injection spectrum of the relativistic par- 
ticles. Following Venter & de Jager (2006) [14], we 
assume that the injection spectrum of the relativistic 
particles Qinj(7,i) obeys a broken power-law 

n U +\ - / ( 2o(*)(7/7br pl for 7 mi „ < 7 < 7b , 

^injlT, t) - | Q oW ( 7 / 7b )-P2 for 7b < 7 < 7max , 

w 

where 7 is the Lorentz factor of the relativistic elec- 
tron and positron. We introduce time independent 
parameters 7 m in, 7b, 7max, Pi and P2 which are the 
minimum, break and the maximum Lorentz factors 
and the power-law indices at the low and high en- 
ergy ranges of the injection spectra, respectively. We 
require that the normalization Qo(t) satisfies the fol- 
lowing equation 

(1 - V )L(t) = /" 7maX Q inj (7, t)7m c c 2 d 7 , (2) 

where m c and c are the mass of an electron (or 
positron) and the speed of light, respectively. 



tial spin-down time of the pulsar. Note that this mag- 
netic field evolution model may be ad hoc, but its be- 
havior is very similar to those adopted in other works 
@, G3- For example, Rees & Gunn (1974) [12] con- 
sider the stretching and winding of the magnetic field 
line inside the uniform PWN and gives B(t) oc t -15 
for n = 3 and B(t) oc £~ 5 / 3 for n — 2.5 for t > r , 
where n is the braking index of the pulsar. Another 
example of de Jager et al. (2009) @ gives B(t) oc t~ 13 
in their calculation of non-relativistic MHD equations. 



C. Evolution of Particle Distribution 

We assume that the distribution of the particles in 
the PWN is isotropic, and then the particle distribu- 
tion function can be easily volume integrated and is 
described by iV(7, t) . The evolution of the particle dis- 
tribution N("f, t) is given by the continuity equation 
in the energy space 

-iV( 7 ,i) + — (7(7, W7,t)) = Qinj (l,t). (4) 

We consider the cooling effects of the relativistic 
particles 7(7, t) including the synchrotron radiation, 
the inverse Compton scattering off the Cosmic Mi- 
crowave Background Radiation (CMB), and the adia- 
batic expansion, i.e., 

7(7.*) = %yn{7,t) + 7ic(7) + 7ad(7.*)- (5) 

Note that the inverse Compton cooling 710(7) does 
not depend on time because we consider that the tar- 
get photon field is only the CMB. As the cooling ef- 
fects, we do not include the inverse Compton scatter- 
ing off the synchrotron radiation field for simplicity 
because it never be a more important cooling process 
than the synchrotron cooling. We include it in the 
calculation of the radiation spectrum. 



B. Evolution of Magnetic Field 

Because the magnetic field lines are stretching and 
winding, it is difficult to model the evolution of the 
magnetic field in the context of the uniform PWN. We 
have to solve the relativistic magnetohydrodynamics 
(MHD) equations to determine the realistic magnetic 
field evolution [6J. For simplicity, we assume that the 
magnetic field evolution are determined in the form of 
the magnetic energy conservation, 

f (i?PWNW) 3 • = J* vHt')dt>. (3) 

In this model, the magnetic field approximately 
evolves as B(t) oc t~ 15 for t > tq, where tq is the ini- 



III. APPLICATION TO THE CRAB NEBULA 

In this section, we apply our model to the Crab Neb- 
ula as the standard calibrator of our model. The Crab 
Nebula is one of the best observed PWN at almost all 
observable wavelengths. 

The Crab Pulsar has the period 3.31 x 10~ 2 s, its 
time derivative 4.21 x 10~ 13 s • s _1 and the braking in- 
dex 2.51. The progenitor supernova is SN1054, which 
means the age of the Crab Nebula t agc « 950yr. Note 
that To ~ 760yr. Here, we adopt that the Crab Neb- 
ula is a sphere of the diameter « 3pc. Combining 
with < agc ps 950yr, the constant expansion velocity is 
vpwn ~ 1500km/s, which is close to the observed ex- 
pansion velocity of the Crab Nebula. 
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FIG. 1: Current spectrum of the Crab Nebula in our model 
and the observational data. The solid line is the total spec- 
trum which is the sum of the synchrotron (dotted line), 
IC/CMB (dashed line), and the SSC (dot-dashed line) 
spectrum, respectively. Used parameters are tabulated in 
Tabled] 



TABLE I: Used parameters to reproduce the current ob- 
served spectrum of the Crab Nebula. 



Adopted Parameter 


Symbol 


Value 


Current Period (s) 


P 


3.31 x 10" 2 


Current Period Derivative (s • s _1 ) 


P 


4.21 x 1CT 13 


Braking Index 


n 


2.51 


Age (yr) 


^agc 


950 


Expansion Velocity (km/s) 




1500 


Fitted Parameter 


Fraction Parameter Eb/(Eb + E e ) 


n 


0.003 


Low Energy Index at Injection 


Pi 


1.5 


High Energy Index at Injection 


1>2 


2.45 


Maximum Energy at Injection 


7max 


7.0 x 10 9 


Break Energy at Injection 


7b 


6.0 x 10 5 


Minimum Energy at Injection 


'"/mill 


1.0 x 10 2 



A. Current Spectrum 

Figure [1] shows the current observed spectrum and 
our calculated one of the Crab Nebula. The adopted 
parameters are shown in Table |TJ As seen in Figure 
III the SSC flux is stronger than the IC/CMB flux at 

7-rays 

The fraction parameter n governs the absolute val- 
ues of the fluxes and the flux ratio of the inverse 
Compton scattering to the synchrotron radiation. The 
KC model obtained u«l from the viewpoint of the 
current dynamical structure of the Crab Nebula, while 
we determine 77 <gc 1 from the viewpoint of the spectral 
evolution. 

In our calculation, the current magnetic field 
strength of the Crab Nebula turns out to be B{t agc ) = 
87^G, which is smaller than B(t agc ) ~ 300/iG used 
by Atoyan & Aharonian [l[. This difference of the 
magnetic field strength can be explained as follows. 



FIG. 2: Spectral evolution of the Crab Nebula. The thin 
solid line is 300yr from the birth. The thick solid, thin 
dotted and the thin dashed lines are lkyr, 3kyr, and lOkyr 
from birth, respectively. Each line represents the total 
spectra which are the sum of the synchrotron, IC/CMB 
and the SSC spectrum. Used parameters are the same as 
in Figure [T] 

Atoyan & Aharonian |l[ adopted B(t agc ) w 300/iG 
from the KC model and adjusted the particle num- 
ber to reproduce the observations. They applied a 
roughly half a spin-down power compared with the 
KC model to reproduce the spectrum, thus another 
half is missing. On the other hand, all the injected 
spin-down power is distributed between the particle 
and the magnetic field energies in our model. 



B. Spectral Evolution 

Figure [2] shows the spectral evolution. As seen in 
Figure [21 the synchrotron flux decreases with time 
and the SSC flux also decreases with time owing to 
decrease of the synchrotron flux, while the IC / CMB 
flux decreases more slowly than the SSC flux. This 
supports the view that old PWNe can be observed as 
7-ray sources which are faint at X-rays. 

The radio/optical observations of the Crab Nebula 
found that the radio/optical flux of the Crab Nebula 
is decreasing with time. The inferred rate of the radio 
flux decrease is — 0.17±0.02%/yr [Hj]. Our model pre- 
dicts the current rate « — 0.16%/yr and this is almost 
consistent with the observation. The inferred rate of 
the optical continuum flux decrease is — 0.55%/yr cal- 
ibrated at 5000 A [l3[ . Our model predicts the current 
rate « — 0.22%/yr, this is by a factor of 2.5 smaller 
than observation. Note that the trend that the de- 
creasing rate increases with frequency is the same as 
the observations. 

To understand the details in Figure [U we show the 
evolution of the particle energy distribution in Fig- 
ure [3) For comparison, the total injected particles till 
lOkyr without the cooling effects is plotted in the dot- 
dashed line in Figure [3] 

In Figure [31 we can divide the evolution of the par- 
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FIG. 3: Evolution of the particle distribution. The thin 
solid line is the distribution at 300yr from the birth. The 
thick solid, thin dotted and the thin dashed lines are those 
at lkyr, 3kyr, and lOkyr, respectively. The dot-dashed line 
is the total injected particles until the age lOkyr. Used 
parameters are the same as in Figure [T] 



ticle distribution in four energy ranges. First, for 
7 > 10 8 , the particle number increases with time. 
Owing to this increase of the high energy particles, 
the IC/CMB flux at lOkyr is larger than that at lkyr 
above 10 TeV. Secondly, for 10 5 < 7 < 10 8 , the evolu- 
tion of the particle distribution is complex. The distri- 
bution is softer than the injection distribution because 
of the synchrotron cooling in this range. Thirdly, for 
10 2 < 7 < 10 5 , the change of the particle distribu- 
tion is small, the difference between the dashed line 
(lOkyr) and dot-dashed line (lOkyr without the cool- 
ing effects) is only a factor of two. This leads to an 
important conclusion that the radio flux decrease is 
mainly because of the decrease of the magnetic field. 
Combining with the observations of the radio flux de- 
crease, our model of the magnetic field evolution can 
be near the truth. Lastly, there exist the particles 
whose energy is lower than 7 m i n . This is because the 
adiabatic cooling is still effective at low energy. 



IV. DISCUSSIONS AND CONCLUSIONS 
A. Discussion 

Our model of the magnetic field evolution is some- 
what ad hoc. The time dependence of the magnetic 
field strength B cx i -1 - 5 for t > To is assumed. This is 
close to other theoretical considerations dill. More- 
over, because our result of the radio flux decrease of 
the Crab Nebula is almost consistent with the obser- 
vation, our model of the magnetic field evolution can 
be near the truth. 

For the injection spectrum of the particle distribu- 
tion, the acceleration of the particles is an unsolved 
problem and we adopt the broken power-law injection. 
It should be noted that one of the important conclu- 
sion in our study that old PWNe can be observed as 7- 



ray sources without a strong X-ray counterpart is not 
affected by the broken power-law assumption. This is 
because the low energy particles do not contribute to 
X-ray and high energy 7-ray emissions. 



B. Conclusions 

In this paper, we built a model of the spectral evo- 
lution of PWNe and applied this model to the Crab 
Nebula as a calibrator of our model. Especially, the 
magnetic field evolution model is unique and can be 
close to the reality. 

The flux decrease of the 7-rays is more moderate 
than radio to X-rays, because the magnetic field de- 
creases rapidly. This result means that old PWNe 
can be observed as 7-ray sources without a strong X- 
ray counterpart. De Jager & Djannati-Atai (2008) [|[ 
also suggested that old PWNe can be observed as TeV 
unidentified sources, because the X-ray emitting par- 
ticles are cooled more rapidly than TeV 7-ray emitting 
particles. This is not the same reason as our conclu- 
sion. 

The current observed spectrum of the Crab Neb- 
ula is reconstructed when the fraction parameter has 
a small value 77 = 0.003. This is consistent with the 
prediction of the magnetization parameter a <C 1 ob- 
tained by Kennel & Coroniti (1984a) @- They ob- 
tained (7 < 1 from the viewpoint of the current dy- 
namical structure of the Crab Nebula, while we deter- 
mine rj <C 1 from the viewpoint of the spectral evolu- 
tion. 

The smaller value of the current magnetic field 
-B(tagc) = 87/iG than inferred value « 300/iG in most 
of other papers is needed to reconstruct the observed 
spectrum of the Crab Nebula. Recent study by Volpi 
et al. (2008) [l|| indicated that the spatially averaged 
magnetic field strength B moa n ~ 100/xG in their rela- 
tivistic MHD simulation. This is close to our value of 
the magnetic field strength. 

Our model can predict the spectral evolution of the 
Crab Nebula, the observed flux decrease of the Crab 
Nebula at radio wavelengths can be explained by our 
magnetic field evolution model. The observed flux de- 
crease rate of the Crab Nebula at optical wavelengths 
is somewhat larger than our model, but the trend 
that the decreasing rate increases with frequency is 
the same as observations. 

Finally, because we can reproduce the current spec- 
trum of the typical PWN with the spectral evolution 
model, we can apply this to other PWNe. We will be 
able to study whether the observed variations of the 
PWNe can be understood by the spectral evolution or 
other effects should be included. 
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